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During the last years it has become clear that global O(N) defects and (7(1) cosmic strings do 
not lead to the pronounced first acoustic peak in the power spectrum of anisotropies of the cosmic 
microwave background which has recently been observed to high accuracy. Inflationary models 
cannot easily accommodate the low second peak indicated by the data. Here we construct causal 
scaling seed models which reproduce the first and second peak. Future, more precise CMB anisotropy 
and polarization experiments will however be able to distinguish them from the ordinary adiabatic 
models. 



PACS: 98.80-k, 98.80Hw, 98.80Cq 
I. INTRODUCTION 

Inflation and topological defects are two classes of 
models to explain the origin of large scale structure in 
the universe. In inflationary models, for fixed cosmolog- 
ical parameters the fluctuation spectrum is determined 
by the initial conditions. In models with topological de- 
fects or other types of seeds, fluctuations in the cosmic 
plasma and in the geometry are continuously induced by 
the gravitational coupling to the seed energy momentum 
tensor. 

Cosmic microwave background (CMB) anisotropies 
provide a excellent link between theoretical predictions 
and observational data. They allow us to distinguish 
between inflationary perturbations and models with de- 
fects by purely linear analysis. On large angular scales, 
both classes of models predict an approximately scale- 
invariant Harrison-Zel'dovich spectrum 0,^]. For in- 
flationary models this can be seen analytically. Scale- 
invariance for defects was discovered numerically ||-^]]; 
simple analytical arguments are given e.g. in ||. 

On smaller angular scales (10'<6><2°), the predic- 
tions of inflation and global O(N) defects are different. 
While inflationary models predict a series of 'acoustic 
peaks', global O(N) defects show a low amplitude broad 
'hump' [^-^]. For local U(l) cosmic strings, the result is 
not so clear. Depending on the detailed modeling of lo- 
cal cosmic strings, the resulting acoustic peaks are quite 
different. The peak can be entirely absent or present 
and even quite substantial, but at an angular harmonic 
£- 400 - 500 

Recent experiments |n| have measured CMB 
anisotropies which are fully compatible with a flat adi- 
abatic inflationary model on the scale of the first peak 
and incompatible with the above mentioned defect mod- 
els. However, the second peak is too low for values of 
the baryon density that are within the constraints in- 
ferred from standard nucleosynthesis. Combining the 



recent BOOMERanG and MAXIMA-I data with infor- 
mations from the distribution of galaxies, a value of 
fl b h 2 = 0.032 ±0.004 was found in M (see also |f| £o|). 



which is incompatible at nearly 3tr with the value f2f,/i = 
0.0189±0.0019 J2^| inferred mainly from measurements of 
primordial deuterium from Ly-alpha absorption systems 
in the continuum emission of 3 high redshift quasars. 

Even if it is fair to say that the possibility of systematic 
errors in all these data sets needs further investigation, 
several, mainly phenomenological, mechanisms have been 
put forward to solve the problem of the low second peak. 
The simplest is clearly to modify standard nucleosynthe- 
sis so that a higher value of the baryon density parameter, 
which leads to a suppression of even peaks becomes ac- 
ceptable p3|-p8t . Another suggestion is to modify one of 
the 'pillars' of the inflationary model, the nearly scale- 
invariant primordial spectrum of fluctuations by adding 
features on it @Q. Also models with Q b h 2 = 0.019 
and with a 'red' tilted spectral index n ~ 0.9, even if 
not preferred with respect to the Vt b h 2 — 0.03 and n = 1 
models, give a reasonable x 2_n t to present data (see, e.g. 
|32| H|). Furthermore, in Ref. H^g^] it has been found 
that a time-varying fine-structure constant can increase 
the compatibility between CMB and BBN data. Finally, 
a combination of inflation with topological defects which 
can contribute to the Sachs- Wolfe plateau and to the first 
peak but not to the second or third peak, has also been 
proposed p8p9| ] as a possible resolution to the problem 
of the low secondary peaks. 

Here we want to investigate whether generic de- 
fect models, the so-called 'causal scaling seeds' mod- 
els, can reproduce the new data. A couple of years 
ago, Neil Turok constructed a model with scaling causal 
seeds which perfectly reproduced the CMB anisotropy 
spectrum of inflationary models j|o| . Other synthesized 
causal seed models with various heights of the acoustic 
peaks are discussed in 41 4^]. Spergel & Zaldarriaga ar- 
gued that causal seeds can nevertheless be distinguished 
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from inflationary models by the induced polarization Jf3[ . 
Our investigations confirm and extend this result. But 
here we shall not only play with some parameters describ- 
ing the model, but we also vary cosmological parameters, 
especially the total curvature which basically determines 
the angular diameter distance and thereby the angular 
scale onto which the peaks in the power spectrum are 
projected. 



II. UNEQUAL TIME CORRELATORS, SEED 
PARAMETERS 



Let us first define the notion of 'causal scaling seeds'. 
Seeds are an inhomogeneously distributed form of energy 
and momentum, which provide a perturbation to the ho- 
mogeneous background fluid. In first order perturbation 
theory they evolve according to the unperturbed (in gen- 
eral non-linear) equations of motion. For simplicity, we 
assume the seeds to be coupled to the cosmic fluid only 
via gravity. A counter example to this are C/(l) cos- 
mic strings. Then the resulting CMB anisotropy power 
spectrum, especially the height of the first acoustic peak, 
depends very sensitively on the details of the coupling of 
string seeds to matter [^T]Q. 

For uncoupled seeds the energy momentum tensor is 
covariantly conserved. To determine power spectra or 
other expectation values which are quadratic in the cos- 
mic perturbations, we just need to know the unequal time 
correlation functions of the seed energy momentum ten- 
sor gj,@, 

(e^(k,» 7 )e; p (k , ,» 7 / )> = M i c lwap (k, v , v ')S(k-k') , 

(i) 

where M is a typical energy scale of the seeds (e.g. the 
symmetry breaking scale for topological defects) which 
determines the overall perturbation amplitude. Seeds 
are causal, if C Ai;y(T p(x, 77, 77') vanishes for |x| > 77 + 77'; 
and they are scaling, if C depends on no other dimen- 
sional parameter than k, 77 and 77'. Using energy momen- 
tum conservation, statistical isotropy and symmetries, 
one can then reduce C^^k, 77, 77') to five functions of 
the variables z 2 — k 2 rji]' and r = rj'/r), which are (a 
consequence of causality) analytic in z 2 EJ. Three of 
these variables describe scalar degrees of freedom, one 
represents vector and one tensor contributions to the 
source correlator C. As in Ref. ||, we parameterize 
the scalar part by the Bardeen potentials of the source, 
e = 4ttGM 2 = 47r(M/M P1 ) 2 , 



(*(k,»7)¥*(k,»/)) = 

($(k,77)<r(k,77')> = 

<*(k,77)$*(k,7/)) = 





(2) 




(3) 


P 3 (z,r) . 


(4) 



The vector and tensor contributions are described by two 
functions S(z,r) and F(z,r) (see Ref. [0 from more 
details). 

Clearly, the parameter space provided by these five 
functions (of two variables) is still enormous and it is 
rather impossible to investigate. For a realistic model, 
the parameter space is even larger due to the radiation- 
matter transition which breaks scale invariance: the seed 
functions can be different in the radiation and in the 
matter era. For global O(N) defects this difference turns 
out not to be very important (less than about 20% |]]) 
it may, however, go to factors of 2 and more for cosmic 
strings Q. 

The topological defect models studied so far, suffer 
from the relatively high amplitude of vector and ten- 
sor perturbations which contribute to the Sachs- Wolfe 
plateau but not to the acoustic peaks. This is the main 
reasons why these models show no significant acoustic 
peaks || . Here, we try to find a causal scaling seed model 
which fits the CMB anisotropy data, hence vector and 
tensor modes have to be suppressed. For simplicity, we 
set £ = F = in this study. In this case, the sum $ + W 
which is due to the anisotropic stresses in the defect en- 
ergy momentum tensor is suppressed by a factor z 2 on 
large scales, z<Cl [Q. In a first attempt we simply set 
* = which implies Pi = P 2 = -P3 = P. 

Another problem of topological defects is 'decoher- 
ence': the coupling of different fc-modes in the defect en- 
ergy momentum tensor, which is due to non-linear evolu- 
tion, 'smears out' distinct features like peaks in the CMB 
anisotropy spectrum into broad humps p7|,p[. To avoid 
this we restrict our study to so called 'perfectly coherent' 
models where the the unequal time correlator P is simply 
the product of the square roots of the two corresponding 
equal time correlators at 77 and 77', 



P(z,r) = \J P(^fz 1 ^, l)P(y/z 2 /r, 1) 



(5) 



This is strictly correct if and only if the time evolution 
of the source is linear. 

In our numerical study described below we investigate 
two families of models. 
Family I 

To enhance the acoustic peak, we use seeds which are 
larger in the radiation era than in the matter era. 



P r {z,l) = 
P m {z,l) 



t 



1 + [bzf 



I 



1 + [bzf 



(6) 



where here the subscripts r and m indicate the radiation 
and matter era respectively. The parameters t and b are 
varied to obtain the best fit and the amplitude e is de- 
termined by the overall normalization. 
Family II 

The second family of models is inspired by Ref. fl40| ], 
which studies spherical exploding shells with p + 3p oc 
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S(r — Arf). To formulate the model we use the source 
functions defined in Ref. |4§| l which determine scalar per- 
turbations of the energy momentum tensor of the seeds, 



e 
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The source functions /« of our models are then given by 
1 sm(Ah]) 



Boo = M 2 f p , 



Qij = M 2 [f p 8 lj - (hkj - ^<Jy)/ir] 



/p + 3/p 



OCT] 



V2 Akf] 



fv = 



k 2 rf>l 2 C 2 



cos(Ckrj) — 



sm(Ckrj) 
Ckri 



with a = {a/a)ri and E = (4 - 2/a)/(3 - 12a). The 
functions f p and f v arc then determined by energy mo- 
mentum conservation fill, 



f P + k 2 f v + -(f p + 3f p ) = 
?7 



V 



fp + ^k 2 U = . 



The function _E is chosen such that the power spectrum 
of f„ is white noise on super horizon scales, a condition 
which is required for purely scalar causal seeds . This 
leads to the Bardeen potentials El] 



r(/p 



A- 2 



3-/„) , 
V 



(8) 
(9) 



Here the seed functions are actually not given as random 
variables but as square-roots of power spectra, and one 
has always to keep in mind that we assume perfect co- 
herence. Of course one can also regard Eqs. (||||) as mere 
definitions with 

AM^fcWA 2 , (10) 

P 2 (z,l) = vk i m 2 /e 2 , (11) 

P 3 (z, 1) = rik^/e 2 = -y/P^z, l)P 2 (z, 1) . (12) 

With a somewhat lengthy calculation one can verify that 
E is chosen such that f n oc const, for z <C 1 and the 
functions Pi(z, 1) are analytic in z 2 — (kr/) 2 . This family 
of models is described by the parameters A and C, which 
have to satisfy < A, C < 1 for causality. Also here one 
can choose different amplitudes for the source functions 
in the radiation and matter era by introduction of the 
additional parameter t ^ 1. 



III. ANGULAR DIAMETER DISTANCE, 
COSMOLOGICAL PARAMETERS 



Seeds generically produce isocurvature perturbations. 
These models, for a flat universe, predict a position of 



the first peak at I ~ 350, which is definitely incompati- 
ble with the recent CMB observations ( po[). How- 
ever, the tight constraints on the flatness of the universe 
obtained from CMB data analysis are based on the as- 
sumption of adiabatic primordial fluctuations. Using this 
loophole, it is possible to construct closed A-dominated 
isocurvature models which have the first acoustic peak in 
the observed position. 

For a given seed-model, the position of the first acous- 
tic peak is determined primarily by the angle subtended 
by the acoustic horizon X ac at decoupling time, % ec . The 
angle under which a given comoving scale A at confor- 
mal time r\dec is seen on the sky is given by 9{X) — 
Vx070 - Vdec), where 

sin(y) if K > 
X (y) = { sinh(y) if K < 
y if K = . 

(K denotes the curvature of 3-space.) 
As the harmonic number I is inversely proportional to the 
angle 0, this yields £ pcak ~ FUg^ where R = 6l?/6 ac . 
The well-known expressions for the conformal time (see 
e.g. Ref. Q) ?7doc and rj are 



K 



^m\/^rad + ^m/(zdcc + 1) 



m -Vdcc = vWk\ 



dz 



o [n m (\+zY + n K {\+z) 2 + n K Y/ 2 



which leads to 



gflat 

ac 



e ac (n m = i,n A = o,n K = o) 

CsVdcc/(m - Vdcc) 

C s V^rad + l/(zdec + 1) . 



where c s = 1/^/3(1 + 3£l/,/4f2 ra d(l + Zdcc)) denotes the 
adiabatic sound speed of the baryon/photon plasma at 
decoupling. We then find 



1 Qrn V^rad + l/(ZfZec + 1) / , „„> 
R = O /LPTr /FT 77 -, _7 ? dec ) ■ (13) 

2 Vl^iflv^rad + "m/^dec + 1) 

Neglecting 51 ra d this reduces to the result of Ref. |52j (the 
factor 1/2 is missing in their formula), 



R = \\lj7TlXino -Vdec)- 



in 



K 



An interesting point is that for Q m — > the quantity 
R depends very sensitively on Thus, we can have 
important shifts in the power spectrum, R ~ 0.6 say, 
with relatively small deviations from flatness (f2 m = 0.3, 
Q A = 0.9, Q K = -0.2). In Ref. ||j the authors have 
shown that the simple prescription t — > Rl reproduces 
the CMB power spectra for curved universes within a 
few percent. On lines of constant R, CMB power spectra 
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are nearly degenerate. In this study we use this simple 
prescription to rescale the flat spectrum. Thereby, we 
make sure that the value of f2 m used in the spectrum 
calculation agrees roughly with the value preferred by our 
best fit value of R and the super- novae constraint [p3[ , 
which can be cast in the form fi m ~ 0.75ft\ — 0.25. f2 m 
determines the time of equal matter and radiation and 
thus influences the early integrated Sachs- Wolfe effect, 
which contributes to the spectrum right in the region of 
the first peak. We therefore get a better approximation 
if we use the correct value for fl m . 



IV. RESULTS 

To analyze family I given by Eqs. (^,0), we have inves- 
tigated a grid of models in (t, b) space with 1 < t < 2 and 
0.1 < b < 1. To make sure that the models are causal, 
we Fourier transform the correlation function into real 
space, cut it at |x| = rj + rj and transform it back. This 
procedure prevents acausal early decay of the correlation 
function; we find that models with b > 1 do not signifi- 
cantly differ from 6=1 after application of this causality 
constraint. 

For each model in our grid we then search the values 
R and the normalization e which minimize \ 2 when com- 
pared with the B98 Q and Maxima Q data. We also 
allow for an overall re-calibration of the B98 data by 20% 
and of the Maxima data by 8%. In Figs. 1 and 2 we show 
the temperature and polarization spectrum for the best 
model (long dashed lines). This model corresponds to 
the best fit parameters t — 2.2, b — 1/9, fl m = 0.35 and 
R = 0.53. It has a value of \ 2 = 38, which, for 22 points 
and 4 parameters (t, 6, R and the normalization), it is 
excluded at more than 99% c.l. if Gaussian statistics is 
assumed. The main disagreement, also for this model, is 
from the high second peak. Assuming flf,h 2 = 0.03 brings 
the model in better agreement with x 2 — 31. However, as 
is clearly visible from Fig. |l|, another main contribution 
to x 2 comes from the last two maxima points. If these 
points are disregarded, the model has a % 2 which is some- 
what lower than the one of a typical ACDM model (short 
dashed line). But it is clearly visible that shifting the 
spectrum does not only move the peak into the correct 
position but it also reduces the width of the peak which 
is already a problem for this model. It is conceivable that 
the introduction of a small amount of decoherence into 
the model might somewhat enlarge the peak width and 
lead to a better fit. Nevertheless, present data already 
does not favor coherent closed isocurvature models over 
the corresponding flat adiabatic models. Furthermore, 
since the model is closed, SIa + £l m ~ 1-2, the secondary 
peaks are at smaller values of £ than in a flat model, which 
makes this model easily distinguishable from a flat model 
with sufficiently accurate measurements as envisaged by 
the Planck satellite J54j . This difference of the inter-peak 
distance which is given only by the values of cosmologi- 



cal parameters like + ^m, is also present in the the 
polarization spectra (see Fig ||). Another important dif- 
ference is that, in general, the signal in the 50 < I < 500 
band is ~ 50% higher for the isocurvature model. CMB 
polarization is produced by Thomson scattering which is 
active only on sub-horizon scales: at fixed £, the relevant 
physical scales are more inside the horizon in the closed 
model and so the contribution to the signal is higher. 



250 




800 



FIG. 1. The CMB temperature anisotropy spectrum 
£(£ + 1)C| T ^ for our best fit model of family I (long dashed, 
blue) and family II (solid, black) is compared with the 
B98 and the Maxima-1 (short dashed, red) data. The 
family I model is a rather good fit to the first peak, 
even if is a closed model (fi ~ 1.2). The family II 
model is flat and is in perfect agreement with the data 
(X 2 = 14/18) even with D. b h 2 = 0.019, as BBN con- 
straints suggest. A standard inflationary spectrum with 
h = 0.65, h 2 Qt — 0.019, flcdm = 0.3, S1a = 1 — is also 
indicated (short dashed). 

A much better fit can be achieved by the models of 
family II. To study these models we have varied 0.3 < 
A,C < 1, and 0.5 < t < 1.5. Our best fit model with 
a value of \ 2 — 14.5 for 22 points and 5 parameters 
(A, C, t, R and the normalization) is in very good agree- 
ment with the data (see fig. [j], solid line), and, up to the 
second peak, is actually quite similar to a model with 
high baryon content. The model shown corresponds to 
the best fit parameters A = 1, C — 0.85, t — 0.8, and 
R = 1. In this model which is flat and causal, the first 
peak in the polarization spectrum is suppressed, as has 
been noted in Ref. 43 (see fig. ||, solid line). 
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FIG. 2. The CMB polarization spectrum C) s for our best 
fit model of family I (long dashed, blue) and family II (solid, 
black) is compared with a standard inflationary spectrum 
with the same parameters as above (short dashed, red). The 
family I model predicts a larger r.m.s. polarization signal in 
the band 50 < I < 500. On the contrary, the lack of interme- 
diate scale polarization at I < 200 in the family II model is 
clearly visible. 

V. CONCLUSIONS 

In this paper we have shown that causal scaling seed 
models for structure formation can reproduce the recent 
CMB anisotropy data |ll|[l6|]. A first very simple closed 
model (family I) can be brought in reasonable agreement 
with all but the last two Maxima- 1 points for a baryon 
density which is not compatible with the nucleosynthesis 
constraint. It is interesting to note that present data al- 
ready slightly disfavors closed isocurvature models since 
they have a peak which is narrower than what is preferred 
by the data. A somewhat more refined model (family II) 
is, for a suitable choice of parameters, in excellent agree- 
ment with all data points in a flat, A-dominated universe. 
The cosmological parameters of our best fit model I are 

n m = n cdm + n b = 0.35, n A = o.ss, h = o.65, h 2 n b = 

0.019 and those of model II are Q m = fl c dm + fib — 
0.4, Q A = 0.6, h = 0.65, h 2 fl b = 0.019. This model is 
preferred by the data with respect to the 'concordance 
model' with fl m = 0.3, Q A = 0.7, h = 0.65, h 2 tt b = 
0.019 and inflationary initial conditions. These models 
can, however, be clearly distinguished from inflationary 
models by future experiments either measuring the sec- 
ondary peaks or the polarization spectrum. The first one, 
a closed model, has smaller inter-peak distances than flat 
inflationary models (see fig. [l]), a definitive lower ampli- 
tude of temperature fluctuations for I > 650 and a greater 



r.m.s. amplitude of polarization for 50 < I < 550; in the 
second one the first peak I ~ 150 in the polarization spec- 
trum is not present (see fig. ^J), which is a consequence 
of causality, and the polarization amplitude is generally 
lower in the band < I < 800. 

To achieve this agreement we have suppressed vector 
and tensor perturbations and have assumed perfectly co- 
herent fluctuations. We believe that it is quite improb- 
able that topological defects from a GUT phase tran- 
sition have such a behavior. Nevertheless, there might 
be some other scale-invariant causal physical mechanism 
(e.g. some spherically symmetric 'neutrino explosions', 
see Ref. p0[) leading to seeds of this or similar type. 
Clearly, we only have a satisfactory model of structure 
formation if also the physical origin of the 'seeds' is clar- 
ified. However, the point of this work was not to find 
"the correct model of large scale structure formation" 
but mainly to investigate, in a phenomenological but at 
the same time physically motivated way, to what extent 
the present values of the cosmological parameters derived 
from accurate CMB data analysis can still be plagued by 
the assumption of the underlying theoretical model. We 
have seen e.g. that flatness, £l m + J1a = 1 is not mainly 
supported by the position of the first peak but by its 
width. Clearly, once secondary peaks are unambiguously 
detected, the inter peak distance will represent another 
direct measure of the total curvature. 

This investigation is rather important especially if 
some of the parameters obtained assuming the standard 
inflationary model are in significant disagreement with 
complementary, more direct observations, as the high 
£lbh 2 value seems to suggest. While present CMB data 
can be regarded as a triumph for a scenario based on pri- 
mordial adiabatic fluctuations, we have presented here 
phenomenological models, based on isocurvature fluctu- 
ations, that also give a good fit to the CMB data. Fortu- 
nately, the concrete models proposed here have peculiar 
characteristics that future CMB experiments will be able 
to detect^]. 
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* A systematic study of the precision with which cosmological 
parameters can be determined by the Planck satellite if the 
model space is enlarged to allow for arbitrary isocurvature 
models without seeds, is presented in Ref. |5q| 
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